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This  is  the  final  report  to  summarize 
some  results  of  FY72  research  task.  The 
investigation  has  attempted  to  provide  a 
basic  understanding  of  plate  impacts. 

Three  categories  of  ballistic  effects  have 
been  considered,  viz.,  perforation  without 
shatter,  perforation  with  shatter  and  no  per¬ 
foration  with  shatter.  Expressions  are 
determined  for  the  ballistic  limit,  residual 
velocity,  shatter  velocity,  and  spray  angle. 
Concepts  from  nonlinear  mechanics  and  shock- 
wave  theory  are  useful  tools  in  this  analys  s, 
and  a  few  nonlinear  differential  equations 
have  been  derived  from  the  projectile's  motion. 
Thus  we  may  speak  of  Bernoulli  nonlinearity  and 
Riccati  nonlinearity.  Our  motivation  is  for 
mathematical  simplicity  with  some  physical 
reality,  and  our  methods  and  results  can  be 
extended  to  more  general  problems  of  penetration 
dynamics . 
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FIGURE 


Figure  1. 


Ballistic  Impact  of  Plate 


INTRODUCTION 


This  is  the  final  report  which  summarizes  some  results  of  my 
investigation  for  the  research  task  of  the  fiscal  year  1972.  The 
investigation  has  been  concerned  with  the  penetration  dynamics  of 
plate  impacts  in  three  ballistic  categories. 

This  research  should  be  of  value  for  consideration  of  weapons 
design  and  armor  protection  from  kinetic  energy  threats. 

A  survey  of  the  literature  [1-4]  shows  that  many  approaches  have 
been  taken  to  provide  solutions  which  characterize  the  ballistic 
impact  with  dominant  features  in  one  way  or  the  other.  But  there 
exists  no  single  approach  which  has  offered  a  convincingly  satis¬ 
factory  solution.  For  this  reason  and  others,  we  seek  to  consider 
the  impact  of  plate  from  the  viewpoint  of  nonlinear  mechanics.  Such 
an  approach  is  especially  suitable  for  the  analytical  treatment  of 
terminal  ballistics  which  is  often  empirical  (and  confidential).  It 
is  not  surprising  that  remarkable  simplification  and  idealization  have 
to  be  made  in  order  to  arrive  at  an  analytical  description  of  the 
problems  considered.  As  mentioned  we  will  consider  normal  impacts 
and  their  ballistic  effects  only.  A  few  modes  of  plate  penetration 

1.  Effects  of  Impact  and  Explosion,  Summary  Technical  Report  of  NDRC, 
Division  2,  Vol.  1,  Chapters  6,  7,  and  p.  279  (Washington,  D.  C., 
1946). 

2.  W.  Herrmann  and  A.  H.  Jones,  Survey  of  Hypervelocity  Impact 
Information,  A.S.R.L.  Report  No.  99-1  (Massachusetts  institute 
of  Technology ,  1961). 

3.  W.  Goldsmith,  Impact ,  pp.  240-248  (Edward  Arnold,  London,  1960). 

4.  N.  Cristcscu,  Dynamic  Plasticity,  Chapter  6  (North-Holi and, 
Amsterdam,  1967). 
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and  perforation  can  be  formulated  as  one-dimensional  nonlinear  prob¬ 
lems.  This  method  we  have  exploited  not  only  for  less  mathematical 
effort  but  also  for  more  physical  support  by  scaling.  Therefore 
our  analytical  results  are  derived  from  concepts  such  as  nonlinear 

resistance  and  shock  fragmentation. 

For  the  categories  of  impact  to  be  considered,  there  are  a  large 

number  of  parameters  which  affect  the  processes  and  results  of  these 
impacts.  Here  we  refer  only  to  the  more  important  factors  such  as  the 
projectile  velocity,  mass  (size,  shape,  density),  tensile  strength, 
plate  thickness,  density,  hardness,  and  shock-compression  properties. 

Upon  impact,  the  projectile  and  plate  interact  with  or  without  perfora¬ 
tion/shatter.  Plate  perforation  may  occur  as  a  result  of  plug  formation, 
hole  enlargement,  bulging  petal,  spallation,  or  fragmentation.  At  the 
same  time  severe  shock  loading  and  unloading  may  cause  large  deformation 
or  complete  disintegration  of  th^  projectile.  Let  us  summarize  these 
situations  as  depicted  (1]  in  Figure  1. 


Plate  Thickness,  in. 

FIGURE  1  -  BALLISTIC  IMPACT  OF  PLATE 


A 


There  we  note  three  lines  of  demarcation:  (a)  the  ballistic  limit 
without  shatter,  separating  impact  categories  I  and  II;  (b)  the  bal¬ 
listic  limit  with  shatter,  separating  categories  III  and  IV;  and  (c) 
the  transition  velocity  dividing  the  no  shatters  and  the  shatters. 
Thus  we  may  consider  the  various  impacts  under  these  headings: 

I.  No  Perforation,  No  Shatter 

II.  Perforation,  No  Shatter 

III.  Perforation,  Shatter 

IV.  No  Perforation,  Shatter. 

In  what  follows,  we  shall  examine  most  of  the  cases  except  for 
category  I  and  petal ing. 

ARMOR  PERFORATION  WITH  PLUG 


Let  us  consider  impact  category  II.  A  simple  model  can  be 
constructed  for  it  on  the  basis  of  nonlinear  resistance  R  *  R(x,x), 
where  x  is  the  penetration  depth  and  x  =  dx/dt  is  the  rate  of 
penetration  with  t  denoting  the  time.  It  is  heuristic  to  consider 
R  =  -  Kxn,  K  and  n  being  appropriate  constants.  Thus  the  equation 
of  motion  of  the  projectile  may  be  written  as 

Mx  +  Kxn  =  0  CD 

where  M  denotes  the  projectile  mass  and  x  its  deceleration.  From 
equation  (1)  the  ballistic  impact  appears  to  be  a  special  kind  of 
nonlinear  vibration.  But  we  will  not  exploit  equation  (1)  in  that 


sense  for  the  projectile's  terminal  trajectory  might  not  be  the  same 
as  the  solution  of  equation  (1).  With  certain  shape  factors  and 
appropriate  energy  balance  accounted  by  constants  K  and  n,  equation  1 
can  satisfy  the  initial  and  boundary  conditions:  (a)  t  =  0,  x  =  0, 
x  =  V  =  striking  velocity;  and  (b)  t  =  T  =  penetration  time,  x  =  P  = 

total  penetration  depth,  x  =  0.  Inserting  these  values  into  the 

first  integral  of  equation  (1),  we  get  the  intial  energy  and  momentum 

as  given  by 


I  MV2  =  KPn+i/(n+l) 

j  i  (n+1) 

MV  =  (2MK/(n+l)r  P 


(2) 

(3) 


Let  P  be  the  plate  thickness.  Then  V  stands  for  the  ballistic  limit 
(V*).  From  Equation  (2)  we  should  have  n  «*  2  for  perforation  by 
energy  absorption  [2]  in  a  volume  proportional  to  P3.  Such  proportion¬ 
ality  is  justifiable  for  pl3te  thickness  of  about  1  caliber;  otherwise 
a  certain  shape  factor  should  be  considered.  From  equation  (3)  a 
momentum- absorption  theory  [2]  requires  o  *  S,  which  implies  also 
the  alternate  form  of  equation  (1)  as 

Mx  ♦  k  x  »  0  (4) 

for  a  different  kind  (k)  of  nonlinear  resistance  (viz.  rate 
sensitive) . 

From  the  similarity  point  of  view,  equation  (1)  can  also  fit 
the  classical  formulations  of  Robins,  Fairbairn,  DoMarre,  and  Krupp 

2.  W.  Hermann  and  A.  H.  Jones,  Survey  of  Hypervelocity  Impact 
Information,  A.S.R.L.  Report  No.  99-1  (Massachusetts  Institute 
of  Technology,  1961). 
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for  the  ballistic  limit  of  armor  plate  [2],  Thus  we  may  put  n  =  0  and 

2 

K^ad  in  equation  (1)  or  (2)  to  derive  the  Robins  formula,  o  being  the 
flow  stress  of  plate  and  d  being  the  projectile  caliber.  For  Fairbaim 
we  should  have  n  =  1  and  K®rd,  t  being  the  shear  strength  of  plate;  for 
DeMarre,  n  =  3/10  and  K°°d^2;  and  for  Krupp,  n  =  1/3  and  K®d^^.  All 
these  different  values  of  n  should  be  attributed  to  the  empirical  nature 
of  the  pertinent  formulas.  As  far  as  terminal  ballistics  and  nonlinear 
mechanics  are  concerned,  our  formulation  of  equations  (1)  and  (4)  is 
analytically  sound.  Slightly  modified  experimental  correlations  are 
up  to  individual  investigators. 

If  the  projectile  has  sufficient  kinetic  energy  to  pass  through 
the  armor  plate,  then  its  residual  velocity  is  given  [5]  by 

l 

Vr  =  M(M  +  in)'1  (V2  -  V*2)2  (5) 

or 

\ 

Vr  =  [M(M  +  m)'1  (V2  -  V*2)]2  (6) 

where  m  is  the  mass  of  the  displaced  plug.  In  equation  (5)  V  should 
be  identified  as  the  flight  velocity  of  the  projectile,  which  differs 
from  the  striking  velocity  due  to  effects  of  friction  and  rotation 
(rifling).  Equation  (6)  neglects  such  difference.  Obviously,  there 
is  no  difference  between  equations  (S)  and  (6)  for  heavy  projectile 
or  thin  plate.  Ideally,  tne  secondary  impact  of  a  double-plate  system 

would  occur  with  almost  the  same  amount  of  momentum  as  the  initial. 

2,  W.  Herrmann  and  A.  H.  Jones,  Survey  of  Hypervelocity  Impact 
Information,  A.S.R.L.  Report  No.  99-1  (Massachusetts  Institute 
of  Technology,  1961) . 

5.  See  Appendix. 
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But  the  impact  energy  is  less.  Threat  from  such  impact  could  be 
indentation,  spalling,  bulging,  or  petaling.  For  a  safe  design  of 
spaced  plates,  Vr  should  be  such  that  the  kinetic  energy  of  the 
secondary  impact  be  small. 

PLATE  PERFORATION  WITH  SPRAY 

This  is  our  Category  III.  See  Figure  1.  For  a  given  projectile 
there  exists  a  shatter  velocity  depending  upon  shock-compression 
properti'  the  impedance  match,  at  and  above  which  it  will  be 

disintegrated  due  to  the  sudden  release  of  strong  shock  pressure 
around  its  free  surface.  Strong  shock  compression  of  solids  requires 
hypervelocity  impact  at  the  level  of  1  km/sec  and  higher  [6] .  We 
shall  give  expression  for  evaluating  the  shatter  velocity  after  we 
consider  the  shock  waves  of  the  explosive  impact. 

Upon  impact  at  hypervelocity,  the  two  colliding  bodies  are 
compressed  by  two  strong  shock  waves  (Sj,  S2)  which  can  be  described 
as  follows: 


(V) 
(8) 
(9) 

6.  Y.  K.  Huang  and  N.  Davids,  J.  Franklin  Inst.  276,  39  (1963). 


°10  +  B  Pl(wi  +  v  -  ui) 

Pi  ■  PhPiCWj  ♦  V) 

Ei  ■  7  pi  ^io'1  -  ’f1) 
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P20W2  "  P2(W2  '  "21 


S2  ?2  ~  P20U2W2 


E2  *  rP2tp20'1  ' 


where  V  is  the  impact  velocity,  W  =  shock  velocity,  u  =  particle 
velocity  behind  the  shock  front,  P  =  density,  p  =  shock  pressure, 

E  =  specific  internal  energy,  and  subscripts  0,  1  and  2  refer  to  the 
initial  value,  striker* and  plate,  respectively.  It  should  be  noted 
that  only  u,  is  a  velocity  relative  to  a  reference  frame  moving  with 

i 

V.  The  interface  between  the  two  colliding  bodies  furnishes  the 
boundary  conditions: 


P2  -  P, 


U2  =  V  -  u^. 


Thus  far  we  have  set  up  only  eight  equations  for  ten  unknowns, 
namely,  p^,  p.^,  Ep  W^,  u^,  p2>  p2>  E2,  W2,  and  u2.  Since  the 
strongly  compressed  material  will  expand  immediately  as  spray,  we 
may  approximate  the  shocked  and  expanded  states  by  a  polytropic 
gas  with  carefully  chosen  index  r  (in  the  limit  r  -*■  3) .  So  two 
additional  equations  are  given  by: 


Ej  *  Pj/Pi  (rx  -  1) 


E2  -  P2/p2  (r2  -  1), 
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Now  we  can  solve  equations  (7) -(16)  simultaneously  for  the  ten 


u2  »  V/U  +  1)  (23) 

W2  =  \  V(r2  +  1)/CC  +  1)  (24) 

P2  =  f  p20(r2  +  1)V2/(C  +  1)2  (25) 

E2  *  jV2/(C  +  l)2  (26) 


Further  results  can  be  sought.  liquations  (15)  and  (16)  imply 
2 

dp/dP  -a  =  rp/p,  a  being  the  sound  velocity.  From  this  and 


yr,r*- 


equations  (17),  (20),  (22),  and  (25)  we  get 


£(5  +  l}~1[r1(r1-H/2]2V 


=  (Wj  +  V)/a]_  =  (rx  +  1)  [2ri(rj  -  1)}  (28) 


(C  +  l)-1[r2(r2  -  1)/2]2V 


=  W2/a^,  =  (r2  +  1)  [2r2(r2  -  1)] 


m  ^ 


with  Hs  denoting  the  shock  Mach  number.  Using  the  Riemann  invariant 
[7],  we  also  get  the  velocity  component  of  lateral  expansion: 


Uj,  =  2a/ (r  -  1) 

the  spray  angle  of  particles  is  given  by  <P  -  tan"1  (u^/u2) , 


Therefore,  we  have 


tan  <f>. 


l  2r  — 

=  5  (— L)2 


,  ,  V  ,  2r2  7 

tan  <f>2  =  -  =  (  - ) 


From  a  or  E  we  can  also  estimate  the  shock  temperature.  In  this 
connection,  the  impact  spray  may  be  thought  of  as  a  result  of 
vaporization  (rather  than  spallation). 

The  foregoing  is  an  analytical  description  for  the  perforation 
of  plate  by  a  small  projectile  shattered  as  a  result  of  shock  loading 


7.  R.  Courant  and  K.  0.  Fredricks,  Supersonic  Flow  and  Shock  Waves 
pp,  87,  90  (Kiley,  New  York,  1948). 
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and  unloading.  It  is  of  interest  to  consider  three  limiting  values 
of  the  index  r.  As  r  4  0,  equations  (18),  (20),  (23),  and  (25)  reduce 
to  these  equations  of  shaped-charge  jet  penetration  [1,2,8).  This 
evidence  seens  to  indicate  the  generality  of  shock-wave  theory  on 
the  one  hand  and  the  compleroentality  of  hydrodynamics  on  the  other. 
Since  two  different  material  models  are  involved  in  the  two  theories, 
we  would  not  seek  the  meaning  of  equations  (15),  (16),  (17),  and  (22) 
with  r  =  0.  A  similar  case  with  r  =  -1  has  been  considered  [7,9] 
for  isentropic  flow  in  gas  dynamics.  Equations  (28)  and  (30)  impose 
another  limit  upon  r.  The  shock  waves  are  said  to  be  hypersonic  with 
Mg  -+•  03  as  r  -*■  1.  Let  r  =  3  and  C  -  1.  From  equations  (32)  and  (33) 
we  get  $  =  60°  which  turns  out  to  be  qualitatively  descriptive  of  the 
spray.  Moreover,  the  pertinent  p,  u,  and  p/pq  are  comparable  with 
those  for  the  products  of  detonation  [10]  (TNT)  based  on  r  -  3.  Thus 
we  consider  r  =  3  appropriate  for  both  dense  gases. 

In  order  to  evaluate  the  shatter  velocity  V  ,  we  must  consider 
some  fundamental  properties  of  shock  waves  in  solids.  It  is  far 
reaching  to  begin  with  the  relation  [11]: 

W  =  A  +  Bu  (34) 

where  A  and  B  are  experimental  constants  with  many  physical  interpreta- 
tions  and  implications.  According  to  Rodean  [12],  A^/2B  is  approxi¬ 
mately  equal  to  the  heat  of  sublimation  of  the  solid.  From  equation 

(18)  the  absolute  particle  velocity  is  given  by  V  -  u^  =  V/ (£  +  1)  and 

?  2 

the  associated  kinetic  energy  is  V  /2(£  +  1)  .  When  the  latter  is 
equal  to  the  heat  of  sublimation,  the  shocked  material  will  expand 
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infinitely  upon  unloading.  Therefore  we  establish  the  relation 

Vg  3  (1  +  Q  A/B  (35) 

1 

T 

where  A  and  B  belong  to  the  projectile  material,  with  £  =  (p2()/pi(P  * 
Current  criterion  for  the  design  of  spaced  plates  is  as  follows. 
At  velocities  where  the  projectile  is  broken  up,  there  is  an  optimum 
shield  thickness  for  which  the  total  penetration  is  a  minimum.  A 
semi-empirical  demonstration  of  such  calculations  is  given  in  Refer¬ 
ence  [13].  At  lower  velocities  with  no  breakup  of  the  projectile, 

the  shield  plate  offers  less  protection  for  potential  threat. 

HYPERVELOCITY  CRATERING 

When  a  target  plate  is  too  thick  to  permit  perforation,  a  crater 
may  form  therein  as  a  result  of  annihilation  of  the  projectile  by 
strong  shock  compression.  Category  IV  of  Figure  1  is  comprised  of 
impacts  like  this.  In  the  literature  there  exist  a  large  number  of 
papers  concerning  this  problem.  See  Reference  [1]  and  Proceedings 
of  Symposia  on  Hypervelocity  Impact  (1958,  1960,  1961,  1963,  1964). 

In  this  report  we  consider  this  problem  again  with  nonlinearities  in 
one  dimension. 

Let  us  denote  the  nonlinear  resistance  by 
R(v,x)  *  a  f(x)v7  ♦  6  g(x)v2 

where  v  is  the  velocity  x  =  dx/dt,  f(x)  and  g(x)  are  two  arbitrary 
functions,  and  a,  8,  and  y  are  all  constants.  Expressing  x  as  v  dv/dx 

1.  Effects  of  Impact  and  Explosion,  Sumnary  Technical  Report  of  NDRC, 
Division  2,  Vol.  1,  Chapters  6.  7,  and  p.  279  (Washington,  D.  C., 
1946). 

13.  l>.  .1.  Stone,  "Designing  for  Meteoroid  Impact,"  ASMF  Paper  67-DE-45 

(March  1%'). 
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with  some  constant  equivalent  mass 
by 

Mv  dv/dx  =  -  a  f(x)vY 


the  equation  of  motion  is  given 
-  6  g(x)v2 


which  may  be  re-arranged  as 


f(x)v’'-1. 

dx  M  M 

Now  equation  (36)  is  of  the  Bernoulli  type.  Let  y  =  v2~Y.  Then 
equation  (36)  can  be  linearised  and  integrated  as 


y  exp  [/  6(2  -  y)g(x)dx/M] 


=  /  (<*/M)(Y  -  2)f(x)  exp  [/  6  (2  ••  v)g(*  )dx/M]dx 

(37) 

It  is  clear  that  we  have  been  motivated  to  formulate  the  idealized 
impact  by  the  well-known  Bernoulli's  nonlinear  differential  equation. 
As  a  special  case,  we  have  the  classical  Poncelet  resistance  [1]  with 
f(x)  =  constant,  g(x)  =  constant,  and  Y  =  0. 

In  the  theory  of  nonlinear  differential  equations,  Riccati's 
equation  is  also  well-known.  It  turns  out  also  reasonable  to  assume 
a  nonlinear  resistance  in  the  form: 

R(v,x)  =  Xv3  -  uvx® 

where  X,  p  and  6  are  all  constants.  Now  the  equation  of  motion 
becomes 

dv  X  y 

—  ♦  -  v2  -  x®  (58) 

dx  M  M 

which  is  the  Riccati  equation.  Let  =  X/M,  C2  *  v/M,  and 

y  =  exp  (Cj/vdx).  Then  equation  (38)  reduces  to 

1.  Effects  of  Impact  and  Lxplosion,  Sumnary  Technical  Report  of  NDRC, 
Division  2,  Vol.  1.  Chapter'  6,  7,  and  p.  279  (Washington,  D.  C., 
1940. 
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0 


d2y 

dx 


2  -ClC2I>'  ’ 


which  is  the  Bessel  equation  with  solution  in  closed  form: 


y  =  x2  [BjI^Chx01)  +  B2I_v(hxcl)] . 
In  equation  (39)  we  have  used: 

v  =  (9  +  2) -1 

l  1 

h  =  2iv(C1C2)2,  i  =  (-1)7 

q  =  —  +  1 

q  2 

4-viri 

Iv(z)  =  e2  J(iz). 

From  equation  (39)  we  get 


(39) 


v  =  Cj’1  d  In  y/dx 

-  (2ClX)-X  [H  2  h  q  +  B2l'v) 

IllVVJ'1!'  (40) 

For  given  initial  and  boundary  values,  we  can  establish  appropriate 
penetration  formulas  from  equation  (40) . 
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CONCLUSION 


This  investigation  has  attempted  to  provide  a  basic  understanding 
of  plate  impacts  from  consideration  of  both  physical  and  mathematical 
nonlinearities.  Cur  physical  concepts  are  nonlinear  resistances  and 
shock  compression  of  solids  -  Our  mathematical  formalism  is  classically 
motivated.  See  equations  (1),  (4),  (36),  and  (38).  Within  this  frame¬ 
work  of  nonlinear  mechanics  and  mathematics,  we  have  examined  ballistic 
effects  of  plate  impact  in  three  categories.  Analytical  evaluation  has 
been  given  of  the  ballistic  limit,  residual  velocity,  shatter  velocity, 
and  spray  angle. 

This  research  is  conceptually  oriented  for  physical  reality  and 
mathematical  simplicity  as  much  as  possible.  Its  methods  and  results 
can  be  extended  to  more  general  problems  of  high-speed  impact. 
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APPENDIX 


Following  Recht  and  Ipson  [14],  we  may  write 

MV  =  (M  +  m)vi  (a) 

i  MV2  -  "  (M  +  m)vi2  =  j  (M  +  m) -1m  MV2  =  Ae  (b) 
with  Ae  denoting  the  energy  loss  due  to  momentum  redistribution.  Let 
AW  =  work  done  for  shearing  the  plug  m  from  the  plate.  Then  energy 
balance  gives 

4  MV2  =  Ae  +  AW  +  j  (M  +  m)Vr2.  (c) 

From  equations  (b)  and  (c),  we  get 

AW*  =  I  (M  +  m)'1  M2V*2  (d) 

for  Vr  =  0  and  V  =  V*  =  the  ballistic  limit.  Assuming  Aw  =  Aw*  and 
substituting  equations  (b)  and  (d)  into  equation  (c),  we  get 
Vr  =  M(M  +  m)"1  (V2  -  V*2)7 
which  is  the  sought  equation  (5) . 

If  we  write  with  AE  *  0 

MV  -  (M  *  m)Vr 

1  MV2  =  AW  +  4  (M  +  m)V  2 

2  2  r 

then  we  obtain 

AW  =  AW*  =  I  MV*2 

and  j 

Vr  =  [M(M  +  m)'1  (V2  -  V*2)]2 

which  is  the  sought  equation  (6) . 

]*•  R-  F.  Recht  and  r.  W.  Ipson,  J.  Appl.  Mech.  30E,  384  (1963). 
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